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Abstract 

We study the difference discrete variational principle in the framework of multi-parameter dif- 
\ ferential approach by regarding the forward difference as an entire geometric object in view of 

■ noncomutative differential geometry. By virtue of this variational principle, we get the difference 

■ discrete Euler-Lagrange equations and canonical ones for the difference discrete versions of the 
classical mechanics and classical field theory. We also explore the difference discrete versions for 
the Euler-Lagrange cohomology and apply them to the symplectic or multisymplectic geometry 
and their preserving properties in both Lagrangian and Hamiltonian formalism. In terms of the 
difference discrete Euler-Lagrange cohomological concepts, we show that the symplectic or mul- 
tisymplectic geometry and their difference discrete structure preserving properties can always be 
established not only in the solution spaces of the discrete Euler-Lagrange/canonical equations 
derived by the difference discrete variational principle but also in the function space in each case 
if and only if the relevant closed Euler-Lagrange cohomological conditions are satisfied. We also 
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^ ■ apply the difference discrete variational principle and cohomological approach directly to the 

■ symplectic and multisymplectic algorithms. 
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1 Introduction 

It is well known that the symplectic structure plays crucially important role in the both 
Lagrangian and Hamiltonian formalism for classical mechanics . On the other hand, 
the multisymplectic structure plays also very important role in the Lagrangian and Hamil- 
tonian formalism for classical field theories [@] IQ ]§] [0 ■ Specially, in the computational 
science, they are extremely important in the symplectic and multisymplectic algorithms 
for the finite dimensional Hamiltonian systems [§0 and infinite-dimensional 010] re- 
spectively. These algorithms are quite powerful and successful in numerical calculations 
of the relevant systems in comparison with other various non-symplectic/multisymplectic 
numerical schemes since the symplectic and multisymplectic schemes preserve the sym- 
plectic structure and multisymplectic structure of the systems respectively. 

Very recently, it has been found |T^[|1^] that there exist what is called the Euler- 



Lagrange cohomology in either classical mechanics or classical field theory and it plays 
very important role for the symplectic or multisymplectic structure preserving property in 
each case. It has also been studied the difference discrete version for classical mechanics 
and field theory mainly in Lagrangian formalism. For this purpose, it has been proposed a 
difference discrete variational principle by regarding the forward (or backward) difference 
as an entire geometric object to deal with variation of the difference discrete classical 
mechanics and field theory [|T^[|lT]. In |]12| and |[T^ [|T^ [|T5| , special investigation has 



been made for the symplectic algorithm as well as the symplectic and multisymplectic 
structure preserving in simple element method respectively from the cohomological point 



of view. In [|I6|, the multi-parameter differential approach has been introduced in order to 
deal with in the same framework the variation of functional and the exterior differential 
calculus in the function space, it has been further studied the Euler-Lagrange cohomology 
and its relation with symplectic and multisymplectic structure preserving properties for 
classical mechanics and field theory in both the Lagrangian and Hamiltonian formalism. 
The cohomological approach has also been applied to what are called Hamiltonian-like 
ODEs and PDEs respectively. 

In this paper, we further study in some details the difference discrete variational 
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principle and apply it to action functional not only in the Lagrangian formalism but 
also in the Hamiltonian formalism. We also study the difference discrete versions for 
the Euler-Lagrange cohomology, symplectic and multisymplectic structures with their 
structure preserving properties in the classical mechanics and field theory. We generalize 
the multi-parameter differential approach for the both variational principle and exterior 
differential calculus in the function space to the difference discrete variational principle 
and exterior differential calculus in the function space on difference discrete base space 
in the Lagrangian and Hamiltonian formalism. It is shown that the difference discrete 
variational principle gives rise to the difference discrete version of the Euler-Lagrange 
equations and that of the canonical equations of motion that preserve the symplectic or 
multisymplectic structures in the Lagrangian and Hamiltonian formalism for the difference 
discrete mechanics and field theory, respectively. It is also shown that the difference 
discrete version of the Euler-Lagrange cohomology in each case is nontrivial and that 
of symplectic and multisymplectic structures are preserved if and only if relevant closed 
Euler-Lagrange conditions are satisfied without making use of the discrete Euler-Lagrange 
equations or the canonical ones in general. Although both the difference discrete Euler- 
Lagrange equations or the canonical ones do satisfy the difference discrete closed Euler- 
Lagrange conditions. Therefore, it is important that these difference discrete version for 
the symplectic and multisymplectic structure-preserving properties hold in the function 
space on the configuration space and its tangent space in the Lagrangian formalism or on 
the phase space in the Hamiltonian formalism with the relevant closed Euler-Lagrange 
conditions in general rather than in the solution space of the difference discrete Euler- 
Lagrange equations or that of the canonical ones only. 

One of the key issues of this paper that is different from the others is the differ- 



ence discrete variational principle is first proposed in [|T^[|rT| to get difference discrete 
Euler-Lagrange equations. As was emphasized in [|l^]||TI|, in view of noncommutative 
differential calculus, the difference is defined as the (discrete) derivative so that it should 
be regarded as an entire geometric object. Furthermore, it can also combine together 
in certain manner as a geometric object to construct the numerical schemes (see section 
5). In the difference discrete variational principle approach, this point of view has been 
carried out. In this paper, this approach is applied not only to the Lagrangian formalism 
but also the Hamiltonian formalism for the both difference discrete mechanics and field 
theory. Together with suitable Leibniz law for the differences, it is also directly applied 
to the derivation of the numerical schemes in symplectic and multisymplectic algorithms. 

Second key issue of this paper is about the difference discrete version of the Euler- 
Lagrange cohomological concepts and content and their role-played in the symplectic and 
multisymplectic structure preserving properties, in each case. As a matter of fact, the 
nontriviality of the difference discrete version for the Euler-Lagrange cohomology plays 
a crucial role and is directly related to the symplectic and multisymplectic structure 
preserving properties. 

In the course of numerical calculation, the "time" t E R is always discretized, say, 
with equal spacing t = At and the space coordinates are also discretized in many cases, 
especially, for the classical field theory. In addition to these computational approach, there 
also exist various discrete physical systems with discrete or difference discrete Lagrangian 
and Hamiltonian functional. It is well known that the differences of functions do not 
obey the ordinary Leibniz law. In order to explore that the difference discrete symplectic 
and multisymplectic structures in these difference discrete systems and their structure- 
preserving properties, some noncommutative differential calculus should be employed. 



3 



Even for the well-established symplectic algorithm. This is the third key point of this 
paper. Recently, the noncommutative differential calculus in regular lattice has been 
employed to deal with the difference discrete phase space for finite dimensional systems 
with separable Hamiltonian |[T^[0[|T9|. Similar noncommutative differential calculus will 



be employed in the present paper. 

Another key point of this paper is the multi-parameter differential approach to the 
difference discrete variational principle and to deal with the exterior differential calculus 
in the function space. This approach provides the same framework for the both DDVP 
and EL cohomological approach. It was employed in ||T6| for the continuous cases. In the 
present paper we will employ the multi-parameter differential approach for the cases on 
the difference discrete version of the base space, i.e. the "time" in the difference discrete 
mechanics and the "spacetime/space" in the difference discrete field theory. 

The plan of this paper is as follows. We first explore, in the framework of the multi- 
parameter differential approach, the difference discrete variational principle in the both 
Lagrangian and Hamiltonian formalism for classical mechanics and classical field theory 
respectively in section 2. It is shown that difference discrete variational principle with 
simply modified Leibniz law for the differences offers the difference discrete version for 
both Euler-Lagrange equations and the canonical equations of motion. In section 3 and 4, 
also in the framework of the multi-parameter differential approach and using the exterior 
differential calculus in the function space in the case of difference discrete base space, we 
deal with such kind of the difference discrete versions of the Euler-Lagrange cohomology 
as well as the symplectic structure preserving and multisymplectic structure preserving 
properties in the both Lagrangian and Hamiltonian formalism for classical mechanics 
and field theory respectively. It is shown that the relevant difference discrete versions of 
the Euler-Lagrange cohomology in each case is nontrivial and it is directly linked with 
the difference discrete symplectic and multisymplectic structure preserving properties. 
We explore in some details the difference discrete variational principle approach and the 
difference discrete Euler-Lagrange cohomological approach to the symplectic and mul- 
tisymplectic algorithms in section 5. It is pointed out that the difference discrete the 
Euler-Lagrange equations, the canonical equations of motion for the classical mechanics 
and field theory present themselves certain symplectic and multisymplectic schemes re- 
spectively. We also show that the Euler midpoint scheme in the symplectic algorithm, the 
midpoint box scheme for a type of PDEs and the midpoint box scheme for the Hamilto- 
nian field theory in the multisymplectic algorithm can be derived by the difference discrete 
variational principle with a suitable difference Leibniz law. And the difference discrete 
Euler-Lagrange cohomology and its relation with the difference discrete symplectic and 
multisymplectic structure preserving properties offer a cohomological scenario to show 
whether numerical schemes are symplectic or multisymplectic. Finally, we end with some 
concluding remarks in section 6. In the appendix, some simple relevant noncommutative 
differential calculus on regular lattice with equal step-length on each direction are given. 
For the sake of self-containing in relevant sections, the content for the continuous case is 
briefly recalled and then the approach is generalized to deal with the difference discrete 
case. 
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2 Variational and difference discrete variational prin- 
ciple in multi-parameter differential approach 

In order to consider certain difference discrete versions of tlie simplectic and multi- 
symplectic structures and their structure-preserving properties in both Lagrangian and 
Hamiltonian formahsm for classical mechanics and field theory, we study the variational 
principle and difference discrete variational principle in both Lagrangian and Hamiltonian 
formalism for classical mechanics and field theory and their difference discrete versions, in 
the framework of the multi-parameter differential approach in this section. We consider 
the cases in classical mechanics in the subsection 2.1, and that in classical field theory in 
the subsection 2.2. 



The difference discrete variational principle approach was first proposed in ||T0| ||TT| with 
vanishing condition at tk = ^oo for the infinitesimal variations of coordinates in the con- 
figuration space, 6q\tk), fc G in the difference discrete classical mechanics and the corre- 
sponding vanishing condition at infinity in 1+1 dimensional or 2 dimensional cases for the 
infinitesimal variations of a set of generic field variables, Su°'^^'^\ a = 1, ■ ■ ■ ,r, {i, j) E ZxZ, 
in the difference discrete classical field theory. As was emphasized, the most important 
point of the approach is regarding the forward difference or its certain combination in 
each difference discrete case as an entire geometric object in the sense of noncommutative 
differential calculus. Of course, if the backward difference is preferred rather than the 
forward one, the framework is almost the same. The framework of the multi-parameter 
differential approach has been employed for the continuous case in |]T^. We review this 



approach and generalize it to deal with the difference discrete variational principle. 



2.1 Variational and difference discrete variational principle in 
continuous and difference discrete classical mechanics 

We begin with recall some content of the multi-parameter differential approach to 
variational principle for classical Lagrangian mechanics and transfer it to the Hamiltonian 
formalism. Then we generalize it to deal with the difference discrete variational principle 
in the both Lagrangian and Hamiltonian formalism for the difference discrete classical 
mechanics. 



2.1.1 Variational principle in multi-parameter differential approach for clas- 
sical mechanics 

Let time t E he the base manifold, M the n-dimensional configuration space on 
t with coordinates q'^{t), {i = 1, • • • TM the tangent bundle of M with coordinates 
{q'' (t) , (t)) , where q^(t) is the time derivative of q^ , F(TM) the function space on TM. 

2.1.1.1 Variational principle in Lagrangian formalism 

The Lagrangian of the systems is denoted by L{q^, q^). For simplicity, we suppose that 
the Lagrangian does not manifestly depend on t. The action functional along a curve q{t) 
in M, with two endpoints a and b, can be constructed by integrating of L along the 
tangent of the curve 

S{q{t)):=J^dtL{q^{t),m)- (1) 



5 



Let us consider the case that at the moment t both and q^{t) variate by an 
infinitesimal increments and the curve becomes a congruence of curves C^J*. The 
infinitesimal variations of and in the congruence can be described as follows 

QKt) ^ qlit) = q\t) + e^5i,q\t), q^{t) qi{t) = q^{t) + eH^q^ (2) 

where e'' , k = 1, ■ ■ ■ ,n, are n free parameters that each of them corresponds one direction 
in the configuration space M, Skq''{t) and Skq^(t) infinitesimal increments of and q^(t) 
at the moment t along the direction k in the congruence of curves 



SkqKt) ■■= ^ 1.^=0 qlit), 6,q^{t) := ^ |,.=o qiit)- (3) 

Here the differentials of g*(t) and qi{t) with respect to in the function space F{TM) 
are manipulated. Namely, the differentials of ql{t) and qi{t) at the moment t in the 
congruence of curves Cea^ can be calculated by: 

dql := ^de' = de%q\ dqi := %de' = deHkq\ (4) 

This framework is called the multi-parameter differential approach. Furthermore, the 
exterior differential calculus in the function space can also be well established in this 
framework. It should be mentioned that in the standard parameter-differential approach 
to the variation calculation is usually to introduce only one free parameter along the curve. 
For the variation calculation it is enough, but it should have more degree of freedom for the 
exterior differential calculation for the functions and functionals. The multi-parameter 
differential setting, in fact, offers the same framework to deal with both variation and 
exterior differential calculation for the functions and functionals. 

In the congruence of curves C^ai the Lagrangian now becomes a family of Lagrangian 
and the same for the action functional: 

sm) - = f,dtLM{t),qm, (5) 

where the upper-index k of e'^ is omitted. 

Hamilton's principle, i.e. the (least) variational principle, seeks the curve Cj' G C^J' 
along which the action S is stationary against all variations of along any directions. 
In the multi-parameter differential approach, similar to @, this can be manipulated by 
taking differentiation with respect to and setting e'^ = afterwards: 

SkSiq{t)):=i,l,=oSMt)) = 0, VA; = l,---,n, (6) 

for all 6ql,{t) = 5kq\t) with 5kq\a) = 5kq\h) = 0. 

It is straightforward to get the differentiation of the action with respect to 

dS.iq.it)) = J^dti,L{ql{t),ql{t))de' = fa dtdq^^ - Iff} + ^dqi ^ ■ (7) 

Therefore, the variation of the action along the direction k is given by 

S.Sit) = J^dt6,q^{§ - If} + §6,q^ \l . (8) 

The last term in the above equation vanishes due to 6kq\a) = 5kq^{h) = 0, hence the 
stationary requirement for S, i.e. the variations of S along any direction should be vanish, 
yields the Euler-Lagrange equations 

aL _ d_dL ^ Q 

dg' dt dq^ 
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2.1.1.2 Variational principle in Hamiltonian formalism 



The action principle can also be carried out on the phase space in the Hamiltonian 
formalism. In order to transfer to the Hamiltonian formalism, we introduce a family of 
conjugate momenta from the family of Lagrangian 

ft. = (10) 
oqi 

and take a Legendre transformation to get the Hamiltonian in the family 

:= H{q\,p,,)=p,,q':-L{q\,qi). (11) 
Now the family of the action functionals can be expressed as 

Se= r dt{pkeq^-H{q\,p,,)} (12) 
Ju 



The variation of the action functional along the direction k can be calculated also in terms 
of differentiation with respect to the parameter e'' and setting e'^ = afterwards 

d f^^ dH dH d 

Thus, the stationary requirement for the action against all variations along any direction, 
i.e. 5kS = 0, V/c = 1, ■ ■ ■ , n together with the fixed endpoint condition lead to the canonical 
equations 

ct = — p= -— (14) 
dpi ' dq^ 

2.1.2 Difference discrete variational principle in multi-parameter differential 
approach for discrete classical mechanics 

Let us now introduce the difference discrete variational principle ||10| |ll[| for the differ- 



ence discrete version of the classical mechanics. Here we generalize the multi-parameter 
differential approach and employ it to deal with the difference discrete variational principle 
for the difference discrete classical mechanics. 

Consider the case that "time" t is difference discretized while the n-dimensional con- 
figuration space Mfc at each moment tk,k G Z, is still continuous and smooth enough. 
Let us assume, without loss generality, that the "time" t & R he discretized as a set of 
nodes and links with equal step-length r = At: 

teR^teT = {{tk,tk+i=tk + r, keZ)}. (15) 

Let A/" and C be the set all nodes and links with index set Ind{N') = Ind{C) = Z, 
M. = Ufcgz the configuration space on T that is still continuous and at least pierce 
wisely smooth enough. At the moment tk, Mk and Ck be the set of nodes and links 
neighboring respectively. For example, includes two links and [t^j^fc+i] 

with endpoints (t^-i, tfc, ^fc+i)- Let the index set of nodes of A4 including t^., Mk = 
Uind{Af)£Tk-^ '^^^ coordinates of are denoted by g*(tfc) = q^^''\i = l,---,n. 
T{Mk) the tangent bundle of in the sense that difference at tk is its base, T*{M^) 
its dual. Let = Uieik^^ union of configurantion spaces M; at ti,l G Ik on 
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A4, TMk = \Jieik ^^io^ of tangent bundles on Mk, F{TMk) and F{TMk) the 

function spaces on each of them respectively, etc.. In the difference discrete variational 
principle, we will use these notions. 

2.1.2.1 Difference discrete variational principle in Lagrangian formalism 

We first study the difference discrete version of the Lagrangian formalism. It is clear 
that the difference discrete Lagrangian written as 

Lz>('=) = Lz>(g*«,gt^«) (16) 

is a functional on F{TAik), since g^^^^) is the forward difference of g**^*^^ at tk defined by 

^^qi{k) _ J^g*^'') = gi*^'') = ijg^e^+i) — g^C')}. (17) 

It is the (discrete) derivative and the base of T(T) in the sense of noncommutative differ- 
ential calculus on a regular lattice with equal step-length r (see also the appendix) 
and the same notation for it as in the continuous case may be employed if it does not 
cause any ambiguity. 

As was emphasized, in what follows the forward difference is viewed as an entire 
geometric object and its dual drt is the base of T*{T) in the sense 

drti^t) = 1. (18) 

It is well known that the (forward) difference as the discrete derivative does not obey 
the Leibniz law but the modified one 

Mf ■ (?) W = Ai/e^) • g('^ + /e^+i) ■ Ai(?W, f,geFM = 0°-. (19) 

On the other hand, however, it is important to note (see the appendix) that in the space 
T*(T) dual to T(T), an exterior differential operator dr exists such that 

dr-.n^T^ dr = 0, (20) 

where the space of Z-forms, / = 0, 1, on T*{T) and dq- does satisfy the Leibniz law: 

d{to A r)(^) = A r^'^) + (-ly^ai^) uj^'') A t^''\ (21) 
The action functional in the continuous case @) now becomes 

SD = T.kezLD{q''^^\qt''^''^), (22) 

where the summation is taken over k & Z. 

We now consider how to calculate the variation of the action functional in this 
case. Since only the "time" is discretized while either the configuration space at each 
moment t^, i.e. at the node fc, or at its neighboring union are still continuous and the 
variational calculation that will be carried out is mainly local, therefore, the difference 
discrete variations may still be manipulated in the framework of the multi-parameter 
differential approach. In addition, as in the continuous case of the classical mechanics, 
the differential and exterior differential calculus in the function space can also be carried 
out in either F{TMk) and F{TMk), etc.. 
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In order to make use of the multi-parameter differential approach the variations of 
qi{k) g^^^ -v^^ith the multi-parameter should be introduced. At the moment tk, we 
have 

qjik) ^ qiik) ^ ^i^^qiik)^ Atg/C^) = Atq''-"^ + e'(5/(Atg*('=)), (23) 

and 

:= |rL^=og/('^ = Siq'^'l SiQu''^'^ := irU=oqtJ'~'^ = Sigt''^"^. (24) 
Then the action functional in ( p2D becomes a family of action functionals 

Sd = J2kez Lde^''^ (25) 

and the variation of the action functional along the direction / 

SiSde = ^5'De|e!=05 (26) 

i.e. 

SiSn. = + ^g^5,g,/W}|,.^o. (27) 

By virtue of the modified Leibniz law (|19D for At = dt, we have 

M^i^^M'''^ = 'i^W^'^ + \i^i^)6,q^i^\ (28) 

Therefore, 

Using the properties (see the appendix) 

Ekez ^tf{tk) = f{tk=+oo) - /(4=-oo), (29) 
and assuming 5iq^^'''^\k±oo = 0, it follows the discrete Euler-Lagrange equations 

It should be mentioned here in general, for the forward difference calculation more gen- 
eral Leibniz law can be adopted and it will lead to more general difference discrete version 
of the Euler-Lagrange equations. We will explore this issue mainly in the Hamiltonian 
formalism in the section 5. 

Let us consider an example. 

Example 2.1. A difference discrete classical mechanics 

Consider a difference discrete version of classical mechanics with following difference 
discrete Lagrangian: 

LB(g^('=), Aig^^'^)) = |(Aig*«)2 - K(g*W). (31) 
The difference discrete variational principle gives the discrete Euler-Lagrange equation 

Ai(Aig*('=-i)) - |V(g*W) = 0, (32) 

i.e. 

9 



This is the difference discrete counterpart of the equation in the continuous case. It has 
correct continuous hmit. 

2.1.2.2 Difference discrete variational principle in Hamiltonian formalism 

Now we consider the difference discrete variational principle on the phase space in the 
difference discrete ( "time" ) Hamiltonian formalism. 

To transfer to the difference discrete Hamiltonian formalism, we first define a family 
of the discrete canonical conjugate momenta 

(fc) _ dLn^'^-^) 

Then a family of the difference discrete Hamiltonian can be introduced through the dis- 
crete Legendre transformation in the family 

i^D.^')=P.('=+^)Ai#)-L^,(^). (35) 

Now the difference discrete version of the action functional in (P^f) becomes a family 
of action functionals: 

Sd ^ Soe = Ekez{p^^'^'^Atq'J'^ - ^De^'^}. (36) 
Then the variation of the action along the direction / in (BH) can be calculated as 



^iS — -^SDe\el^=0- (37) 

For the differential of Sd^ in the above equation, we have 

dSn. = j:.,z{dpJ'^'\A,gi^''^ - f^) - {A,p^ + '-^)de^ + A,{v^ dc^^^'^))} . 

(38) 

Here, the modified Leibniz law ([T9| ) has been used. 

Now, the difference discrete variational principle gives rise to the difference discrete 
version of the canonical equations of motion ([T^) 

^^^f"" = f*S^> ^^P^''' = (39) 

In fact, the first set of equations above can directly be derived from the Legendge 
transformation (p5[ ) and the second set can be gotten from the Legendge transformation 
([35| ) and the Euler-Lagrange equations (|30|) . This indicates that the difference discrete 
variational principle approach to the difference discrete version of classical mechanics is 
self-consistent. 

Let us consider the example 2.1 in the discrete Hamiltonian formalism. 
Example 2.2. Hamiltonian formalism for the example 2.1. 
First, the difference discrete conjugate momentum is introduced 

{k) ^ Mp'"-^' _ A,o^(^-i) r4nl 
P' a(Vt<?"('=-i)) ~ *^ • ^^^1 

The Hamiltonian is introduced through the discrete Legendre transformation 

i/^W = pXk+^)/\^qm - L^(fe) = ip/^+i)' + V(g^W). (41) 
And a pair of difference discrete canonical equations read now 

In fact, the time difference discrete derivative of pi'-'^-* can also be derived from the differ- 
ence discrete Lagrangian and the discrete Euler-Lagrange equation (B^). 
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2.2 Variational and difference discrete variational principle in 
continuous and difference discrete classical field theory 

We now study the difference discrete variational principle in Lagrangian and Hamilto- 
nian formalism for the difference discrete classical field theory. For the sake of simplicity, 
let us consider the 1+1-d and 2-d cases in discrete classical field theory for a set of generic 
fields m", a = 1, ■ ■ ■ ,r. We first recall the multi-parameter differential approach to the 
variation of functional in Lagrangian formalism and deal with the Hamiltonian for- 
malism, then generalize it to the difference discrete variational principle for the difference 
discrete classical field theory in both Lagrangian and Hamiltonian formalism. 

2.2.1 Variational principle in multi-parameter differential approach for clas- 
sical field theory 

For the sake of simplicity, let X*^^'"~^) be an n-dimensional Minkowskian space as base 
manifold with coordinates x'^, (yU = 0, ■ ■ ■ , n — 1), M the configuration space on 
with a set of generic fields ^"(x), {a = l,---,r), TM the tangent bundle of M with 
coordinates where = F{TM) the function space on TM etc. We also 

assume these fields to be free of constraints. In fact, the approach here can easily be 
applied to other cases. 

2.2.1.1 Variational principle in Lagrangian formalism 

The Lagrangian of the fields now is a functional of the set of generic fields under 
consideration: 

m°) = / d"-ia;/:(M"(x, t),u'^{^, t)), = u"(x, t), etc., (43) 

and the action is given by 

S{u") = J dtL{u", m") = / d"a;/:(u"(a;), u'^ix)), (44) 

where C{u°'{x),u'^{x)) is the Lagrangian density. 

In order to apply Hamilton's principle we first consider how to define the variation 
of the action functional S{u°') in a manner analog to the case of classical mechanics. In 
order to achieve this purpose, let us suppose that both u^i^x) and u'^{x) variate by an 
infinitesimal increments such that at a spacetime point of x the infinitesimal variations 
of u°' and -u" can be described as follows 

M^(x) = m"(x) + e^5/3M°(a;), u'^^{x) = u'^{x) + e^5/3U°(x), (45) 

where e'^, P = 1, ■ ■ ■ , r are free parameters, each of which corresponds one direction in the 
configuration space M, and 

6f3U^{x) := £j Ip^o <(a;), ^Kei^) ■= A l^''=o <e(a;), (46) 

the infinitesimal increments of u°'{x) and u'^{x) along the direction f3 at the spacetime 
point X respectively. 

On the other hand, similar to the case of classical mechanics, the differential calculus of 
u'^{x) and u'^{x) in the function space F{TM) can also be manipulated in the framework of 
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the multi-parameter differential approach. Furthermore, the exterior differential calculus 
in this frameworlf can afso be weff estabhshed. 

Thus the (exterior) differentials of u°'{x) and u'^{x) in the function space F{TM) at 
the spacetime point x can be defined as: 

du'^ix) := ^^def^ = def^Spu'', du%{x) := ^^de^ = de^Spu'^^. (47) 

Now, the Lagrangian also becomes a family of Lagrangian functionals 

L,«,«^) = /c^"-ia;£«(x,t),M^,(x,t)), (48) 

and the action S{u°') becomes a family of functionals as well 

S^S,^S{uf). (49) 

Then the variation of the action along the direction (3 can be manipulated as the derivative 
of yS'g with respect to and setting = afterwards. Namely, 

SpS := ^ Se. (50) 

Manipulating the variation of the action functional in this manner and integrating by 
parts, it follows that 

SpS = /d"a;{(|| - d,i^J)Spu^ + a^(^<5^<)} . (51) 

Assuming 6pu'^ |±oo= 0, and requiring 6pS = along all directions according to Hamilton's 
principle, then the Euler-Lagrange equations follow: 

B - 9.i-§^J - 0. (52) 

2.2.1.2 Variational principle in Hamiltonian formalism 

In order to use the multi-parameter differential approach for the variational principle 
in Hamiltonian formalism for the classical field theory, we first define a family of a set of 
"momenta" that are canonically conjugate to the family of field variables 

^^eix) = (53) 

and take a Legendre transformation to get the Hamiltonian density in the family 

K«, TTae, V„0 = 7r„,(a:)^/°(x) - A«, iif, V„0, (54) 

where Va = a = 1, ■ ■ ■ ,n — 1. A family of the Hamiltonian then is given by 

H,{t)^ Jd^-^xUeix), (55) 

with the Legendre transformation 

H,{t) = / d^-^XT:Ux)u''^{x) - L,{t). (56) 

The action S{u°') becomes a family of functionals as well 

S^S, = j (i"xKe(^)«^(^) - V„<)}. (57) 
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Then the variation of the action along the direction (3 can be manipulated as the derivative 
of with respect to and setting = afterwards as was shown in Namely, 

(58) 

Thus, the canonical equations of motion follow form the stationary requirement of the 
action principle 

^"(^) = ^^(^) = -B + (59) 

2.2.2 Difference discrete variational principle in multi-parameter differential 
approach for discrete classical field theory 

For the sake of simplicity, we consider the case of 1+1-dimensional spacetime or 2- 
dimensional space. It is straightforward to generalize for higher dimensional case. 

Let X*-^'^) or X^"^^ with suitable signature of the metrics be the base manifold, = X 
a regular lattice with 2-directions x^, (/x = 1,2) on X*-^'^-* or X^'^\ Af the all nodes on 
that are coordinated by X(^ij), E Z x Z with index set Ind{M), Md := M(ij) the 
pierce of configuration space with a set of generic field variables M"(a;(jj)) = g Md 

at the node TM(^ij) the tangent bundle of M(j j) with the set of field variables and 

their differences {u"^''^\'u'^^''^^) G T(M(ij)), F{TM(iJ)) the function space on TM^ij^, etc.. 

For a given node with coordinates let A/(jj) be the set of nodes neighboring to 

with index set X(jj) = Ind{N)[ij), = Uind{j\r)£i(^i ^^-^ ^ set of nodes that is 

related to by the differences, A4d '■= -Mx^.j-, = iJind{Af)€X(i -^J^ ^^e union of the 
pierces of configuration space on X(^ijy 

It is known that the forward differences along each direction in F{TA4 j,) are defined 
by A^m(^'^') = 

Aim(*'^') = ^^{u'-'+^'^^ - «(*'■'■)), AaM^*'-'') = )^(m^''^'+^^ - u^''^'^). (60) 

They are the bases of T{X) and the upper-indexes reflect the corresponding coordinates 
of nodes on X. And their dual rfx^ = d^x^ are the bases of T*{X) 

dxxf'iA,) = 5'i. (61) 

As in the previous subsection, the (forward) differences as the discrete derivatives do 
not obey the Leibniz law but the modified one (|1^) along each direction. While in the 
space T*{X) dual to T(A'), an exterior differential calculus can be introduced (see the 
appendix) such that there exists an operator dx with the following properties 

dx : d\ = 0, (62) 

where is the space of all /-forms in T*{X) and dx does satisfy the Leibniz law: 

dx{uj A r)(^'J') = dxuj^'^^^ A t^''^^ + {-ly^ai^^u^^'^) A dxT^''^\ (63) 

It is important to note that although the base manifold is discretized either the con- 
figuration space at each node or its neighboring union is still continuous. In addition, 
the variational calculation that will be carried out is mainly local. Therefore, similar 
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to the case of difference discrete classical mechanics, in what follows the difference dis- 
crete variations will be manipulated on the framework of the multi-parameter differential 
approach. 

In addition, as was pointed out in the difference discrete version of classical mechanics, 
the differential and the exterior differential calculus in the function space can also be 
carried out in either F(TM(jj)) or F{TA4x^i j^,) ■~ UaTgj -^(^■^)7 ^tc. on the framework 
of the multi-parameter differential approach. 

2.2.2.1 Difference discrete variational principle in Lagrangian formalism 

For the difference discrete version of the classical field theory, the difference discrete 
Lagrangian denoted as 

= (64) 
is a functional in F{TA4x^i j^)- The action functional is given by 

Sd = Ei^J)eZxzCD{u''^'''\uf'^^). (65) 

Taking the variation of by the multi-parameter differential approach, the variation 
along the direction (3 is given by 

SfsSo = ^SDe\e0=O- (66) 

For the sake of simplicity, in what follows we omit the multi-parameters in the course 
of calculation. Thus, we have 

Employing the modified Leibniz law ([T9|) for the forward difference, we have 

Assuming that 5/3m"*^'^'''*'s vanish at infinity, it follows the discrete Euler-Lagrange equa- 
tions 

Let us consider an example to show that the difference discrete variational principle 
gives right result. 

Example 2.3. A discrete classical field theory with difference discrete Lagrangian 
Consider the following discrete classical field theory with difference discrete Lagrangian: 

Cd{u''^''^\u'^^''^^) = i(AX^''^'^)^ - V'Im"^''^'^)- (68) 

The difference discrete variational principle gives the discrete Euler-Lagrange equation 
([67| ) as follows 

Ai(Aiu"(^'i'J')) + A2(A2?i"(*'-''"^)) - ^l^(n°(^'^')) = 0, (69) 
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I.e. 

(70) 

This is also what is wanted for the difference discrete counterpart of the relevant Euler- 
Lagrange equation in the continuous limit. 

2.2.2.2 Difference discrete variational principle in Hamiltonian formalism 

Let X^^'^^ be the base space. We first define a set of the discrete canonical conjugate 
momenta on the tangent space of the set of nodes neighboring to the node x{i,j), i.e. on 

The difference discrete Hamiltonian is introduced through the discrete Legendre transfor- 
mation 

Hd^^''^ = 7r„(*+i'^')AtM"(^'^') - Cd^'''\ (72) 
The action functional (^) now is given by 

Sd = E(,,)e^xz(vr.(^+^'^')A,w"(^'^-) - Hd^^^'^)- (73) 

Taking the variation of 5*^ by the multi-parameter differential approach and using the 
modified Leibniz law, the variation along the direction (3 is given by 

S,Sn = E(,,)ezxz{5,vr«(^+i'^-)(A,«"(^'^-) - £^0^) - {A,nJ^'^^ + ^ 

(74) 

+ A,(vr«(^'^-)5,n"(^'^-)) - A.(^(£$;^5,«"(^-^-))}. 

The stationary requirement for difference discrete version of action functional against all 
variations along any direction gives rise to the difference discrete version of the canonical 
equations of motion: 

(75) 

_ _^!1d L A (_£J1d ^ 



It should be noted that the first set of equations above may directly follow from the 
difference discrete Legendre transformation ([72|), while the second set from the transfor- 
mation (|72D and the difference discrete Euler-Lagrange equations (^). As in the case of 
difference discrete version of the classical mechanics, this also indicates that the difference 
discrete variational principle approach is self-consistent. 

Let us consider the example 2.3 in the discrete Hamiltonian formalism. 

Example 2.4: Difference discrete Hamiltonian formalism for the example 2.3. 

First, a set of the difference discrete conjugate momenta are introduced 

Then the Hamiltonian is introduced through the discrete Legendre transformation 
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And a pair sets of difference discrete canonical equations read now 

A,«-(^'^-) = 7r"(^+i'^-), A,7r„(^'^-) = -fjS^. (78) 

In fact, the time discrete derivative of tTq,^*'-'-' follows from the difference discrete Lagrangian 
and the discrete Euler-Lagrange equations (|69D. 



3 Euler-Lagrange cohomology, symplectic structure 
preserving property in continuous and discrete clas- 
sical mechanics 

Now we are ready to study certain difference discrete version of the Euler-Lagrange 
cohomology, its relations to the simplectic and multisymplectic structures as well as their 
preserving properties in the classical mechanics and field theory in both Lagrangian and 
Hamiltonian formalism. In this and the next section, we deal with the case of classical 
mechanics and field theory respectively. We first recall some content on the issues in the 
continuous case for the mechanics in subsection 3.1. Then we consider these issues in 
difference discrete classical mechanics in the Lagrangian and Hamiltonian formalism in 
subsection 3.2 and 3.3 respectively. 

3.1 Euler-Lagrange cohomology, symplectic structure preserv- 
ing property in continuous classical mechanics 

By virtue of the multi-parameter differential approach to exterior derivatives in the 
function space, the Euler-Lagrange cohomology and symplectic structure preserving prop- 



erty in Lagrangian formalism for the classical mechanics has been studied in ||T6[ . We recall 
the content in 3.1.1 and deal with the the relevant issues in the Hamiltonian formalism 
in 3.1.2. 



3.1.1 Lagrangian formalism 

It is obvious but important to note that in manipulating the variation of the action 
functionalfor the classical mechanics in Lagrangian formalism in the subsection 2.1, the 
variation of the action is very closely linked with the differential of with respect to the 
free parameters e'^ in the congruence of curves, i.e. dS^- 

Furthermore, the exterior differential calculus in the framework of the multi-parameter 
space can be well established. Namely, in the free multi-parameter space, the standard 
exterior differential calculus can be introduced. And this multi-parameter exterior differ- 
ential calculus can also be employed to deal with the exterior differential calculus in the 
function spaces on M and TM, i.e. FM and F{TM). 

Therefore, the integrand in (0) for dS^ and the boundary term may be views as 1- 
forms. And the differential of the family of Lagrangian functions := L{ql{t), ql{t)) with 
respect to is given by 

dL, |..=o= m - If + i{§dql}. (79) 
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Let us define a family of the Euler-Lagrange 1-forms and that of the canonical 1-forms 
e on T*M, 

^efea:={ff-|||W, (80) 
Ol. := ^dql, (81) 



we have 

(ni ni\ I J^fni ni\ J , 

dt 

Furthermore, owing to the nilpotency of d with respect to e'^ in the cotangent space 
of the congruence of curves on F{T*M), 

d^LM,Qi) 1.^=0=0, 

it follows that 



dE{q\q') + iu;L = 0, (83) 



A 

dt 

where ul is the symplectic 2-form in the Lagrangian formalism defined by 



= d9L = ^dq^ A dq^ + ^dq^ A dq\ (84) 

And it does not change if the canonical 1-form transforms as 

e^eL' = eL + d(3{q\q'), (85) 

where (3{q^,q^) is an arbitrary function of {q\q^)- 

We have established the important and significant issues on the Euler-Lagrange co- 
homology and symplectic structure preserving law in the Lagrangian formalism. Before 
enumerating them and exploring their significance, let us investigate the relevant issues 
in the Hamiltonian formalism. 

3.1.2 Hamiltonian mechanics 

The above cohomological and other issues in the Lagrangian formalism of the classical 
mechanics can also be well established in the Hamiltonian formalism. Again the multi- 
parameter differential approach will be employed. 

In order to transfer to the phase space of the Hamiltonian formalism, in the subsection 
2.1, we introduce a family of conjugate momenta (jT^) from the family of Lagrangian L^, 
i.e. 



dqi 

and taken a Legendre transformation to get a set of the Hamiltonian functions in the 
family in (|ll|) , i.e. 

H,:=H{q\,p,:)=PkA'^-L{q\Al). 

We have also calculated the variation of action functional along the direction k in ([T3|) 
where the differentiation of a family of action functionals have been taken as follows 

dS, = dtM - lt)dp,, - (If + P^e)dql + iMM)}- (86) 

Let us introduce a pair sets of canonical Euler-Lagrange 1-form families in the Hamil- 
tonian formalism as follows: 

E,M^V,.) = -{P,e + l^)dq\, E,M\.Pi.) = (Qi ' I^Jdpj,- (87) 
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And a family of the canonical 1-forms 

One = Piedqi (88) 

Thus, the equation (|86|) can be expressed as 

dS, = \ll dt{E^, + E,, + pHe}. (89) 

This equation can be viewed as the equation for 1-forms. 

Furthermore, due to the nilpotency of d in the multi-parameter differential approach, 
it is straightforward from 

= ^25-^1^, = (90) 



to get 

dt 

Here ujh is the symplectic 2-form in the Hamiltonian formalism 



d{Ep + E^) + jrUH = 0. (91) 



duH = dOn = dpi A dq\ (92) 

We may introduce a family of zj = {pj , qj), where pj = {pi^, ■ ■ ■ , Pne), qj = {q^e: ' y ^"e 
defined by 

z,{t) :=z{t)+e'6iz{t), (93) 

where 

6iz{t) = ^ |,.=o ze{t) (94) 

infinitesimal variation of z{t) along the direction / in the configuration space. 
Then the Euler-Lagrange 1-forms in (^) become 

E,,{z,, 4) = dzf{V,,H, - Jz,), (95) 

where J is a 2r;, x 2n symplectic matrix. 

Taking the exterior differential of E^e and setting e'^ = afterwards, it follows that 

dE, + \dz'^ A Jdz = 0. (96) 

This is equivalent to the equation (pi]). 

3.1.3 Remarks on Euler-Lagrange cohomology and symplectic structure pre- 
serving property in classical mechanics 

It is easy to check that under the Legendre transformation (^) all relevant issues in the 
Lagrangian and Hamiltonian formalism are in one-to-one correspondence. For instance, 
we have the following equivalent relations: 

Lagrangian formalism Hamiltonian formalism 

L{q\ q') ^ pi(f - H{pi, q^) 

E{q\ct) ^ Epip„q^) + E,ip„q^) 

Ol ^ en (97) 

dL = E + pL ^ dipif -H) = Ep + Eg + pH 

LOl ^ UJh 

dE + PL = ^ diEp + Eg) + Ph = 
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Therefore, these important issues on the Euler-Lagrange cohomology and the sym- 
plectic structure preserving property in two formahsms are corresponding to each other. 
Let us enumerate these issues in the Hamiltonian formahsm. 

First, the null canonical Euler-Lagrange 1-forms with e' = give rise to a pair sets of 
the canonical equations of motion (0). In terms of z, the null canonical Euler-Lagrange 
1-form ( ^5|) with = gives rise to the canonical equations in z{t): 

z = J-^V,H. (98) 

Secondly, the above null forms are the special case of the coboundary canonical Euler- 
Lagrange 1-forms, say, 

E{z,z) = da{z,z), (99) 

where a{z^ z) is an arbitrary function of (z, z) on the phase space. 

Thirdly, from either the expression of d{piq^ — H) or the definitions of the canonical 
Euler-Lagrange 1-forms it is easy to see that the canonical Euler-Lagrange 1-forms are 
not exact in general. Therefore, there exists a nontrivial Euler-Lagrange cohomology in 
the Hamiltonian formalism for classical mechanics defined as 

-f^CM^ ={Closed Euler-Lagrange forms}/ {Exact Euler-Lagrange forms}. 



Fourthly, from the equation (|91|) or ( ^6]) it follows a theorem on the symplectic struc- 
ture preserving law in the Hamiltonian mechanics. 
Theorem 3.1: 

On the phase space of the Hamiltonian mechanics, the symplectic structure uih pre- 
serving law 

i^H = (100) 



holds if and only if the canonical Euler-Lagrange forms are closed: 

dE{z,z) = 0, I.e. d{Ep{q\p,)+E,{q\p,))=0. (101) 



This means that the symplectic conservation law holds not only in the solution space of 
the equations of motion as shown in the standard approach [|T| but also in the function 
space with the closed Euler-Lagrange condition ||T0|1 [pJ]] . 

Finally, it should be mentioned that the Euler-Lagrange cohomological scenario in the 
Hamiltonian formalism may be performed in two slightly different processes. Namely, it 
may either start from the exterior derivative of the action functional or begin with the 
canonical equations. In the second process, the families of the canonical Euler-Lagrange 
1-forms ( p7|) and (|95D may be introduced directly referring to the corresponding canonical 
equations and (p8| ) but releasing first all canonical variables from the solution space of 
the canonical equations ([I^ ) and (|98D respectively. Then by taking the exterior differential 
of the families of the canonical Euler-Lagrange 1-forms (^) and and setting the free 
parameters being vanish afterwards it also follows the theorem 3.1 as was shown in (pG]). 

3.2 Discrete Euler-Lagrange cohomology and symplectic struc- 
ture preserving property in difference discrete classical me- 
chanics 
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3.2.1 Difference discrete Lagrangian formalism 

Let us first consider tlie discrete Euler-Lagrange coliomology, its relation to tlie sim- 
plectic structure and its preserving property for tlie difference discrete classical mechanics 
in the Lagrangian formalism. 



The difference discrete Lagrangian at the moment tk is written as in (16) on F{TAik)- 
That is 



i{k) i(k) 



1 1 



where qt^^^ is the forward difference of g**^'^-' defined in (|T]_ 

Taking the exterior differential d of L/)'^'^^ in the function space F{TAik) as has been 
done in the subsection 2.1 in the framework of multi-parameter differential approach, it 
follows 

Using the modified Leibniz law (^) with respect to the forward difference = dt defined 

in ( p!?!) and introducing the discrete Euler-Lagrange 1-form as well as the discrete canonical 

(fc) 

:= _ A*(ag#Sy))}^^^^'^ (102) 

(^DL^'' = M^^)dq^'''^ (103) 

we have 

(ILd'^^^ = Ed^^'^ + Ai^^(*^\ (104) 
Due to the nilpotency of d on T*{Mk), d'^Lo^^'' = 0, we get 

dEn^''^ + AtUoL^"^ = 0, (105) 



1-form Oj;)]^ 



where ujdl"^^ is a discrete symplectic 2-form on T*(A^ 



k 



(106) 

3.2.2 Difference discrete Hamiltonian mechanics 

Let us now study the symplectic structure preserving property in Hamiltonian formal- 
ism for the difference discrete mechanics when the time t is discretized in the manner 
of last subsection. We start with what has been constructed for the difference discrete 
mechanics in Lagrangian formalism. 

Let us first define the discrete canonical conjugate momentum 

Pi^^'' = afAtg'C^-i))- (10''') 

Now the difference discrete Hamiltonian can be introduced through the discrete Legendre 
transformation 

Hn^'^ = pi('=+i)Aig^W - L^^^^). (108) 

A set of canonical equations for the time difference discrete derivative of pf'^ follow from 
the difference discrete Euler-Lagrange equations (|30|) and the above discrete Legendre 
transformation (|108|) 

Am^'^ = (109) 
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From the discrete Legendre transformation, another set of canonical equations for the 
time difference discrete derivative of g**^^-* follow 

= (110) 

In terms of z^^^ , a pair of the canonical equations P^ ) become 

A,;^^ = J-^\/,Hd^''\z'^''^). (Ill) 

We now consider the relevant cohomological issues in the difference discrete version 
of the Hamiltonian mechanics. This time we start from introducing the discrete Euler- 
Lagrange 1-forms. To this end, all time discrete canonical variables should 
be released from the solution space of the difference discrete canonical equations (ff^) or 



(|111|) . This can easily be realized by means of the multi-parameter differential approach. 



Thus a pair of discrete Euler-Lagrange 1-forms can be introduced: 

fll2l 



or in term of z^''^ 



Edz^^'KzC'^) = rf;2«'^(JAtzW - V,Hd^''\z^''^)). (113) 



By taking the exterior differential of the discrete Euler-Lagrange 1-forms ( p.l2| ) and 
( |113| ), it is straightforward to prove the following formula 

dEn.^'\z^''^) + AtcunH^''^ = 0, (114) 

where udh^'^^ is the difference discrete version of the symplectic 2-form at the moment tk 
given by 

cu^^W = irf^W^A Jd2(^). (115) 

We can also start from the exterior differential of L^,^''^ = pj^^'^+^^A^g**^*^'' — Hd'^'^\ 
introduce the discrete Euler-Lagrange 1-forms and the canonical 1-form then take the 
second exterior differential to get above equations. 

3.2.3 Remarks on discrete Euler-Lagrange cohomology and symplectic struc- 
ture preserving property in classical mechanics 

As in the case of continuous classical mechanics, it is easy to check that under the 
difference discrete version of the Legendre transformation ( |108| ) all relevant issues in the 
difference discrete Lagrangian and difference discrete Hamiltonian formalism are in one- 
to-one correspondence. In fact, we have the following equivalent relations: 

Discrete Lagrangian formalism Discrete Hamiltonian formalism 

AT (k) _ p (k) j_ dn (fc) ^ rl ( ^^''^'^'^ Ai(k) ZTik)\ _ pik) , p{k) , d n{k) 

dEo^"^ + iouDL^'^ = ^ d{EDp^'^+EDg^''^) + iu;DH^^^ = 

(116) 
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Therefore, for enumerating the important issues on the difference discrete Euler- 
Lagrange cohomology and the difference discrete symplectic structure preserving property, 
we may also work with one formahsm, the corresponding issues in another formahsm can 
easily be established. Let us also enumerate these issues in the difference discrete Hamil- 
tonian formalism as follows. 

First, the null discrete canonical Euler-Lagrange forms give rise to the canonical equa- 
tions and they are the special case of the coboundary discrete canonical Euler-Lagrange 
forms. 

Secondly, since the first terms in the definitions ( [L12| ) and ( |113| ) are not exact in general 



so that the discrete canonical Euler-Lagrange forms are not always exact. Therefore, there 
exists a nontrivial difference discrete version of the Euler-Lagrange cohomology in discrete 
Hamiltonian mechanics: 

HDCM'-={dosed discrete Euler-Lagrange forms}/{Exact discrete Euler-Lagrange forms}. 
Thirdly, from the equation (|114|) it follows straightforwardly the following theorem 



for the necessary and sufficient condition of the difference discrete symplectic structure 
preserving law. 
Theorem 3.2: 

The difference discrete symplectic structure preserving equation 

AtUDH^''^ = 0, I.e. u;dh^^+^^ = ujdh^''^ (117) 
holds if and only if the discrete Euler-Lagrange forms are closed: 

dEDz^^\z^^^) = 0. (118) 



Fourthly, the difference discrete symplectic structure preserving law holds in the func- 
tion space associated with the difference discrete version of the closed Euler-Lagrange 
condition in general rather than in the solution space of the canonical equations only. 

Finally, all these issues can be reached by taking the exterior differential of pi^''^^^ q'^^''^ — 
Hd^'^^ first, then introducing the discrete Euler-Lagrange 1-forms and discrete canonical 
1-form. The theorem 3.2 follows from (i^(pj('^+^)g*('^) —Hd^'^^) = 0. In fact, in the difference 
discrete version of the Lagrangian formalism, it has started and progressed in this manner. 



4 Euler-Lagrange cohomology, multisymplectic struc- 
ture preserving property in continuous and differ- 
ence discrete classical field theory 

We now study the discrete Euler-Lagrange cohomology and its relation to the discrete 
multisymplectic structure in discrete classical field theory. In order to self-contained, 
we first recall the the content on Euler-Lagrange cohomology and multisymplectic struc- 
ture preserving property in continuous classical field theory |TI[] as well as the multi- 



parameter differential approach to the exterior differential in the function spacein the 
subsection 4.1 for the both Lagrangian and Hamiltonian formalism. For the sake of sim- 
plicity, for the difference discrete cases, let us consider the l-|-l-d and 2-d cases in discrete 
classical field theory for the a set of generic fields without constraints in the Lagrangian 
and Hamiltonian formalism respectively in the subsection 4.2. 
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4.1 Euler-Lagrange cohomology, multisymplectic structure pre- 
serving property in continuous classical field theory 

We recall some content on the Euler-Lagrange cohomology and its relation to the 
multisymplectic structure preserving property in both Lagrangian and Hamiltonian for- 
malism for continuous classical field theory by means of the multi-parameter differential 
approach. 



4.1.1 Classical field theory in Lagrangian formalism 

In the multi-parameter differential approach to the variation of the action of the clas- 
sical field theory in subsection 2.2.1, the variation along the direction (3 of the action has 
been given by the equation (^0]), i.e. 



d_ 



^13^ '■— TTJ^ |e/5=0 ^e- 



It ia shown that the differentiation of the action functional with respect to the free pa- 
rameters is given by 

Note that the integrand in the above equation is an equation for 1-forms: 

= (H - M§il-))d^l + M^du-). (119) 
Let us define a family of the Euler-Lagrange 1-forms 

u%) := - d,{§^))du-, (120) 
and n sets of 1-forms that each set corresponds to a family of canonical 1-forms 



■= (121) 



Then the equation ( |119[ ) becomes 



dC, = E,{u^^,u%) + d^e^,. (122) 

Furthermore, due to the nilpotency of d with respect to e^, taking the second exterior 
differential of C^{u",u^^) and setting = afterwards d'^ C^{u" ,u^^) 1^/3=0= 0; it follows 
that 

dE{u^,u^^) + £,UL^ = Q, (123) 
where ujl^ are n symplectic structures defined by 

= = ^^du^ A du" + ^§;pdu^, A duf" . (124) 
And they do not change if the set of n canonical 1-forms transform as 

^ e^'' = r + rf/3(u°, u'^), (125) 
where /?(«", m^) is an arbitrary function of (m",m^). 
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4.1.2 Classical field theory in Hamiltonian formalism 

All steps that have been progressed in the Lagrangian field theory can also be well 
progressed very similarly in the Hamiltonian formalism. In order to do so, we first define 
a set of "momenta" that are canonically conjugate to the field variables 

M^) = 

and take a Legendre transformation (Q) to get the Hamiltonian density 

n{u", Tia, Vau'') = 7i^{x)u''{x) - C{u'' , m", V^m"), a = 1, ■ ■ ■ , n - 1. 

Let us consider the action S'(M")becoming a family of functionals as follows 

S^S, = /rf"x{7r,,(x)M^(x) - 7^,«,7r,„ VaO}. (126) 

The differential of S'^ can be calculated in the multi-parameter differential approach and 
the integrand's differential dC^ reads 

(127) 

+Vt{n^du^) - Vaio^^^du"). 

Introducing a pair sets of the canonical Euler-Lagrange 1-forms in the Hamiltonian 
formalism 

E^ = dn^{u-{x)-^^}, E. = rf«-{-7r„(x)-^ + V.^(^^}, (128) 
and n-canonical 1-forms 

00 = TT^du", = Q^^)du", (129) 

then we have 

dC.U^o = E^ + E^ + Vt^o - Va^^ (130) 

Due to the nilpotency of d, (i^£e|e/3=o = 0, it is straightforward to get the following 
formula: 

d{E^ + E^) + Vtc^o - V,c^" = 0, (131) 
where uj^ and uj'^ a set of n symplectic 2-forms 



- rfvr, A du^, uj^ = d(^-^i^) A du^. (132) 



4.1.3 Remarks on Euler-Lagrange cohomology and multisymplectic structure 
preserving for classical field theory 

It is easy to find that, similar to the case of classical mechanics, under the Legendre 
transformation (Q) all relevant issues in the Lagrangian and Hamiltonian formalism are 
in one-to-one correspondent equivalent. In fact, we have also the following equivalent 
relations: 



Lagrangian formalism Hamiltonian formalism 

E{u^,u^") ^ Eu + E^ 



Ol^ ^ eu^eH (133) 



dC = E + O^Ol^ ^ d(7r„M" -n) = Eu + E^ + O^Oh^ 
dE + df.ujL^' = ^ d{Eu + K) + d^oju" = 
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Therefore, for enumerating the important issues on the Euler-Lagrange cohomology 
and the multisymplectic structure preserving property in classical field theory, we may 
work with one formalism, then the corresponding issues in another formalism are indi- 
cated automatically. Let us this time enumerate and verify these relevant issues in the 
Lagrangian formalism. 

First, if the Euler-Lagrange 1-form in ( |120|) is null with = 0, i.e. 

1.^=0=0, (134) 

it gives rise to the Euler-Lagrange equations (|5^). 

Secondly, E^{u",u'^^ = is a special case of the coboundary Euler-Lagrange 1-forms 

E,«,M-J =rf«,«,M-J, (135) 

where a^{u'^,u'^^) a family of arbitrary functions of (M",M^e)- Although they are cohomo- 
logically trivial but it can already be seen that in the Euler-Lagrange 1-forms, (m",m^) 
are already not in the solution space of the Euler-Lagrange equations only rather they 
are in the function space with corresponding closed Euler-Lagrange condition (see below) 
in general. 

Thirdly, if the Lagrangian density C in (|48|) changes to C by adding certain term 



£(m", m°) ^ C{u^, M^) = £(«-, M") + V(m"), (136) 

where V(m") is an arbitrary function of u", the equation ( |119| ) changes to 

dC[l,=,= E'{u^,u'^) + ^e^, (137) 

where E^u"", m^) differs from E{u°', m^) by changing C to C in the expressions, while a set 
of n canonical 1-forms 9^ have not been changed because V(m") does not depend on m^. 
In fact, the Euler-Lagrange equations have been changed by adding a potential-like term 
that does not depend on u'^. This means that even if by adding a coboundary term, the 
Euler-Lagrange equations do change and the set of canonical forms may still be the same 
as before. Furthermore, the canonical transformations or the multisymplectic mappings 
that preserve the Euler-Lagrange equations, correspondingly the canonical equations in 
the Hamiltonian formalism, and Euler-Lagrange 1-forms as well as the multi-symplectic 
structures form invariant. This will lead to the issues on the generating functions, its 
relation to the Euler-Lagrange cohomology and so on. We will explore these issues else- 
where. 

Fourthly, from the equation (|119|) it is easy to see that i?e(M",M" ) in general are not 



cohomologically trivial because the families of canonical 1-forms are not trivial. Therefore, 
There exists a nontrivial Euler-Lagrange cohomology in the classical field theory for the 
set of generic fields u'^{x): 

HcFT'-={ closed Euler-Lagrange forms}/{ exact Euler-Lagrange forms}. 



Furthermore, From the equations ( |123| ) and ( p.31| ) in the Hamiltonian formalism it 
follows straightforwardly an important theorem in the classical Lagrangian field theory. 
Theorem 4-1: 

For a given Lagrangian field theory, there exists a set of n symplectic structures u^'^ 
and the multisymplectic preserving property, i.e. the conservation or divergence free law 
of the multisymplectic structures 

£:^L^ = (138) 
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holds if and only if the relevant Euler- Lagrange 1-form is closed 



dE{u'',u'^) = 0. 



(139) 



It is interesting to see that if we introduce a new 2-form 



fi(M°,M^) = rfE(M-,M-). 



(140) 



Then Q may be viewed as a [/(l)-hke curvature 2-form while the Euler-Lagrange 1- 
form the ?7(l)-hke connection 1-form. Therefore, the closed Euler-Lagrange condition is 
nothing but the flat connection condition. On the other hand, if for some reason that the 
multisymplectic conservation law is broken then the broken pattern may be described by 
the curvature 2-form Q. There is a similar issue in the finite dimensional case as well. 

It is also important to notice that the multisymplectic structure preserving property 
is directly linked with the closed Euler-Lagrange condition. And although the null Euler- 
Lagrange 1-form, the coboundary Euler-Lagrange 1-forms satisfy the Euler-Lagrange con- 
dition, it does not mean that the closed Euler-Lagrange 1-forms can always be exact as 
was pointed out above. In addition, u"(x)'s in the Euler-Lagrange condition are not in the 
solution space of the Euler-Lagrange equations only in general. Therefore, the multisym- 
plectic structure preserving property, i.e. the conservation law of the set of n symplectic 
2-forms u^, holds not only in the solution space of the field equations but also in the 
function space with the closed Euler-Lagrange condition in general. 

Analog to the case of classical mechanics, there is also an another slightly different 
way to deal with the issues on the Euler-Lagrange cohomology and the multisymplectic 
structure preserving property. Namely, either taking the exterior differential of the La- 
grangian first, or directly starting from the Euler-Lagrange equations ([5^). By means of 
the multi-parameter differential approach, it is easy to release all field variables form the 
solution space of the Euler-Lagrange equations to the function space F{TM) and intro- 
duce a family of the Euler-Lagrange 1-forms ( |120| ) associated with the equations (^21). 
Then by taking the exterior differential of the Euler-Lagrange 1-forms and setting e° = 0, 
it is straightforwardly to re-derive the theorem on the necessary and sufficient condition 
for the multisymplectic structure preserving law. 

4.2 Discrete Euler-Lagrange cohomology and multisymplectic 
structure preserving property in discrete field theory 

For the sake of simplicity, we consider the cases of 1+1 or 2 dimensional base manifold. 
Let X^'^ or X"^ with suitable signature of the metrics be the base manifold, a regular 
lattice with 2-directions x^, (/x = 1,2) on X^'^ or X^, the configuration space with 
g and so forth as before. 

4.2.1 Discrete classical field theory in Lagrangian formalism 

The difference discrete Lagrangian for a set of the generic fields w", a = 1, ■ ■ ■ , r, is a 



functional in F{T{M.x^^ ■^)) 




(141) 
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Taking exterior differential d e T*{A4x^iji) of C^^^'^^ and making use of tlie modified 
Leibniz law (0), in the framework of multi-parameter differential approach, we get 

dCD^'^^^ = EDiu"^'''\ AX^^'^')) + A,9dl^^''''\ (142) 

where E^^'^'^^ are the discrete Euler-Lagrange 1-forms defined by 

and Odl^^^'''^ are two canonical 1-forms: 



It is easy to see that there exist two symplectic 2-forms on T*{Aix^^ j-,)'- 

cu^^Mm) = rf^^^MM)^ /i = l,2. (145) 

The equation d'^Co^^'-'^ = 0, on T*{A4x^, j-,) leads to the discrete multisymplectic structure 
preserving property, i.e. the conservation law or the divergence free equation of cu^^^'^^: 

dED{u'^^''^\ A^m"(^'^)) + A^ujdl"^''^'^ = 0. (146) 

4.2.2 Discrete classical field theory in Hamiltonian formalism 

Let us now study the difference discrete classical field theory in Hamiltonian formalism 
in the case of the spacetime/space x^,^ = 1,2, are discretized in the manner of last 
subsection while the configuration space at each node of the relevant lattice with 
coordinates x{i,j), E Z x Z and the ones on its neighboring are still continuous. 

In order to transfer to the difference discrete version of the Hamiltonian formalism, 
we first define a set of the discrete canonical conjugate momenta on the tangent space of 
the set of nodes neighboring to the node x{i,j), i.e. on T(A'(jj)): 



° a(Aiw"(»-i.j))- \^^' ) 



The time difference discrete derivative of vr"'*^*'-'^ follow from the discrete Euler-Lagrange 
equation (|32| ) 

Now the difference discrete Hamiltonian is introduced through the discrete Legendre trans- 
formation 

Ud^^'^^ = 7r„(^+^'-')AtM"(^'^') - Cd^''^'>. (149) 
The difference discrete version of action functional is given by 

SdH = jCd^'''^ = Ei^,,){rcJ'''''''>^tU"'^'''^ - Hd^''^^ ■} (150) 

The differential of each term can be calculated in the multi-parameter differential approach 
dCr. (*'^')| n r. - dTT (»+i.i)i'A,7/"(*'J') - MLd^1£L) _ (A,7r ^^'^'^ 4- 
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(151) 



Now we introduce the difference discrete version of the canonical Euler-Lagrange 1- 
forms 



and the difference discrete version of two canonical 1-forms 



(152) 



Thus the equation (|151| ) becomes 

dCDe^'''\,=, = EdJ'''^ + EdJ'''^ + Mdh''^'^'^ - AJdh"^'''^. (154) 

Now due to the nilpotency of d, d'^CDe^^'''^\ei^=o = O5 it follows the equation for difference 
discrete multisymplectic structure preserving property: 

dEoJ'''^ + dEnJ'''^ + AtUDH^'^'''^ - A.ljdh'^'''^ = 0, (155) 

where lodh^^'^'^^ and ludj^^^^'^^ are two symplectic 2-forms: 



4.2.3 Remarks on discrete Euler-Lagrange cohomology and multisymplectic 
structure preserving property in discrete field theory 

It is clear that similar to the case for the classical field theory we have also the fol- 
lowing equivalent relations between the difference discrete versions of the Lagrangian and 
Hamiltonian formalism for the difference discrete classical field theory: 

Discrete Lagrangian formalism Discrete Hamiltonian formalism 

En'^''^^ ^ Ed.'^'''^ + EdJ''''^ 
dCo^''^^ = El^'^^^ + A^Old^^'^^^ ^ d{^^u^ - 7^)(^'^) = {Edu + + A^0,5i/^)(*'^) 
dEo^'''^ + A^,ujLD^^'^^^ = ^ d{E^ + E,)g'^') + A^ujhd''^'''^ = 

(157) 

Let us enumerate and verify relevant important issues on the discrete Euler-Lagrange 
cohomology and difference discrete multisymplectic structure preserving property in the 
difference discrete Lagrangian formalism. 

First, the null discrete Euler-Lagrange 1-form corresponds to the discrete Euler-Lagrange 
equations and it is a special case of coboundary discrete Euler-Lagrange 1-forms 

^^(^J) = daD^''^\ (158) 

where a^)*^*'-'^ an arbitrary function on F{T*Aix!^,j-))- 

Secondly, although they satisfy the discrete Euler-Lagrange condition, it does not 
mean that all closed discrete Euler-Lagrange 1-forms are exact. As a matter of fact, from 
the equation ( |142| ) it is easy to see that the Euler-Lagrange 1-forms are not exact in 
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general since the two canonical 1-forms 6j:)^^^'^\ {jj, = 1,2) are not trivial. Therefore, for 
a given difference discrete field theory, there exists a nontrivial difference discrete version 
of the Euler-Lagrange cohomology: 

-f^DCFT:= {closed Euler-Lagrange forms}/{exact Euler-Lagrange forms}. 

Thirdly, from the equations (|146|) and ( |155|) it follows the theorem on the necessary 
and sufficient condition for the difference discrete multisymplectic structure preserving 
law: 

Theorem 4-2: 

The difference discrete multisymplectic structure preserving law 

^,uji''^ = (159) 

holds if and only if the discrete Euler-Lagrange 1-form satisfies the discrete Euler-Lagrange 
condition, i.e. it is closed: 

dED^''^^ = 0. (160) 



In addition, this also indicates that the variables ■u^'^'^'^^'s etc. in the cohomology are 
still in the function space rather than the ones in the solution space only. Consequently, 
this means that the difference discrete multisymplectic structure preserving law holds in 
the function space with the closed discrete Euler-Lagrange condition in general rather 
than in the solution space only. 

Finally, it should be mentioned that all these issues can be straightforward to generalize 
to higher dimensional cases of spacetime or space X" as base manifold. 



5 Difference discrete variational principle and dis- 
crete Euler-Lagrange cohomological approach to 
symplectic and multisymplectic algorithms 

It is worthwhile to show that the difference discrete variational principle and the coho- 
mological scenario described in the previous sections for the difference discrete mechanics 
as well as for the discrete classical field theory in the both Lagrangian and Hamilto- 
nian formalism can be directly applied to the numerical schemes in the symplectic and 
multisymplectic algorithms respectively. 

As a matter of fact, the difference discrete versions of the Euler-Lagrange equations 
and the canonical equations of motion in the both Lagrangian and Hamiltonian formalism 
offer themselves the numerical schemes in symplectic and multisymplectic algorithms for 
the difference discrete version of the classical mechanics and the difference discrete version 
of the classical field theory respectively. 



In addition, as has been pointed out in that the Euler-Lagrange-like cohomological 



approach may also be applied to the so-called Hamiltonian-like ODEs and PDEs respec- 
tively. Therefore, for the difference discrete versions of these Hamiltonian-like ODEs and 
PDEs the difference discrete variational principle and the difference discrete version of 
the cohomological approach may also be available. 

We investigate the relevant issues for the symplectic algorithm first in the subsection 
5.1, and the relevant issues in the multisymplectic algorithm first in the subsection 5.2. 
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5.1 Difference discrete variational principle and discrete coho- 
mological approach to symplectic algorithm 

As in the case of both Lagrangian and Hamihonian formahsm for classical mechanics, 
the above-established difference discrete variational principle and the discrete cohomolog- 
ical scenario in the difference discrete Lagrangian and Hamiltonian formalism for classical 
mechanics should be directly applied to the difference discrete versions of ODE's with La- 
grangian and/or Hamiltonian and to the numerical schemes for the symplectic algorithm. 
In this section we study this topic. We first consider how to apply the difference discrete 
variational principle to the numerical schemes in symplectic algorithm in the subsection 
5.1.1. Then we show the difference discrete version of the Euler-Lagrange coho mo logical 
scenario offers how to justify a scheme is symplectic in subsection 5.1.2. 



5.1.1 Difference discrete variational principle and symplectic algorithm 

Let us first note that both the difference discrete versions of the Euler-Lagrange equa- 
tions (|30D in the Lagrangian formalism and the difference discrete versions of the canonical 



equations (^2[) or (|111|) in the Hamiltonian mechanics are derived from a relevant differ- 
ence discrete variational principle in the framework of the multiparameter differential 
approach. These difference discrete equations may in fact offer certain numerical schemes 
and are automatically symplectic. 

In the difference discrete Lagrangian formalism case, the discrete Euler-Lagrange equa- 
tions (|30f) read 



Their offering a numerical scheme that is symplectic can be seen manifestly from the 
difference discrete Euler-Lagrange equations (|3^) in the example 2.1: 



oq 

i.e. 

oq 

For given initial values of g*('=^=o)^ g*(fc=i)^ ^j^e equations give rise to g'C^^^) and so on. For 
the Lagrangian ( |3TD this is in fact a simplest scheme that preserves, in the sense of time 
difference discrete, the symplectic structure. This has been proved in the subsection 3.1. 

For the case of difference discrete Hamiltonian mechanics , the difference discrete 
canonical equations in the (^2]) or (|1 1 1|) read 



and 

They also offer a set of numerical schemes in difference discrete Hamiltonian formalism 
and preserve the symplectic structure in the sense of time difference discrete. This has 
been proved in the subsection 3.1 as well. 
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In fact, in the multi-parameter differential approach to difference discrete variational 
principle, it is easy to see why the difference discrete variational principle offers the nu- 
merical schemes that are automatically symplectic. It has been shown that the variations 
of the difference discrete action functionals are calculated by taking the differential of the 
difference discrete action functionals and setting e'^ = afterwards, i.e. dLo^l^k^Q in the 
Lagrangian formalism and d{piAtq'^ — H)D^\^k^Q in the Hamiltonian formalism. As was 
shown in the subsections 2.1 and 3.1, the stationary requirement of difference discrete 
variational principle leads to either difference discrete Euler-Lagrange equations or differ- 
ence discrete canonical equations of motion. Furthermore, the second exterior differentials 
of the difference discrete action functionals lead to the symplectic structure preserving 
law. 

Let us consider the Euler mid-poind scheme. We find that there exists a difference 
discrete action functional for the scheme. 

The difference discrete Lagrangian of the scheme can be given by 

Lmdpt^^^ = Atg^W), ^^^^+5) := ^g^C^+i) + g^W). (161) 

Note that at the moment tk q^^'^^^^ and Ajg*'^'^^ are the coordinates and the tangents for 
the scheme. 

Then the discrete canonical momenta conjugated to the coordinates can be defined as 

„ (k+i) _ dLrr^dptW ._ 1/ (fc+1) (fc)x 

The discrete Legendre transformation is given by 

= p,('=+i)Aig*W - L,ndpt^'\ (163) 

The difference discrete action functional of the scheme is given by 

SmM = E.ez{pr'/'^A,g^W - H{q^('^^/'\pt'^'^}. (164) 

Let us take the variation of this difference discrete action functional in the framework 
of the multi-parameter differential approach, i.e. the variation along the direction / is 
manipulated by 

^iSmdpt = '^Smdpte\e'-=Q- (165) 

Eventually, the differential of each term under the summation of Smdpt is given by 
:= rfp,(^'+i/2)(A^^.W _ ^H^q^ik+l/2)^pXk+m^) 

(166) 

-{AtPi^''^ + ^i7(g^('=+V2),p.(fc+i/2)))rfg*(fc+i/2) + A^(p.(fc)c;g^(fc)). 

Note that here the following generalized modified Leibniz law (see also the appendix) 
has been used: 

At(/W . ^W) := i(/('=+i) . ^(^+1) - /W . ^W) 

(167) 

= -{Atf^''^ ■ ^(^=+1/2) + fik+im . Ai^(*=)). 

The first equality in the above equation is the definition for the forward difference, while 
the second can be easily proved. 
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It is now clear that the stationary requirement of the difference discrete variational 
principle leads to the Euler mid-point scheme as follows: 

(168) 

Thus we have shown that the Euler mid-point scheme is a difference discrete variational 
scheme with corresponding (dependent) difference discrete variables and suitable Leibniz 
law for the difference. 



5.1.2 Discrete cohomological approach to symplectic algorithm 

We now consider how to apply the difference discrete version of the Euler-Lagrange 
cohomological approach to the numerical schemes in symplectic algorithm. 

In the standard approach, it is commonly accustomed to regarding a numerical scheme 
as a (time-discrete) mapping. In order to justify whether a given numerical scheme is 
symplectic, the standard approach is to see whether this mapping is symplectic preserving 
and the verification is always carried out in the solution space of the scheme In the 

difference discrete version of the Euler-Lagrange cohomological approach, however, instead 
of working on the solution space, it is working on the function space with the relevant 
cohomological issues. Analog to the case for the difference discrete classical mechanics, 
there are two slightly different ways to apply the cohomological approach. Namely, one 
is based upon the difference discrete variational principle for the schemes and taking 
second (exterior) differential of the action functional to get the necessary and sufficient 
condition for the symplectic structure preserving property of the scheme. Another is to 
release the scheme away from the solution space, even if the solution space does exist, 
and to introduce some suitable difference discrete Euler-Lagrange 1-forms associated with 
the scheme such that the null difference discrete Euler-Lagrange 1-forms give rise to the 
scheme. Then by taking the exterior derivative of the difference discrete Euler-Lagrange 
1-forms to see whether follows a time-discrete symplectic structure preserving law. 

We will consider some examples to show how the cohomological scenario works. Let us 
first retain to the mid-point scheme ( |168| ) and follow-up by the cohomological approach. 

It is clear that by introducing the discrete Euler-Lagrange 1-forms 



(169) 



EpC^) = -(A^p.W + |iJ(g^('=+i/2),p.(fc+i/2)))rfg^(fc+i/2)^ 
and the difference discrete canonical 1-form for the scheme: 

Orndpt^'^ =P^^''^dt('^\ (170) 

the equation ( |166| ) becomes 

dL^,,J'^ = E,^'^ + E,^'^ + AtO^mdpt^'y (171) 

By taking exterior differential of the above equation ( |1 71| ) , due to d^L^dpt = 0, it follows 
that 

d{E^ + ^p)^^) + Atiu^d^/'^^ = 0, (172) 
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where oJmdpt"'^'^ = dpi^''^ A dq^'^^\ Since the null Euler-Lagrange 1-forms of ( |169|) give rise 
to the scheme and automatically satisfy d{Eq + Ep)^^'> = 0, while the latter leads to the 
symplectic conservation law from (|172|) : 

A , , e^) ■— , C^+i) , , (>')\ — n 

^t^mdpt ■— —y^mdpt — ^mdpt ) — U. 

Therefore, the midpoint scheme is symplectic. 

On the other hand, we can also start directly from the the mid-point scheme in terms 
of z^''^: 

AjzW = J-iV,//(i(z(^+i) + (173) 

This time we release the scheme form the solution space first, even if it does exist. This 
can be treated by the multi-parameter differential approach. For simplicity, we suppose 
this has been done already. Then we introduce a difference discrete Euler-Lagrange 1-form 
for the scheme 

E.mdpt^"^ = ldiz(>'+'^ + JA,^« - V.H{l{z('+') + zW))} (174) 

such that the null discrete Euler-Lagrange form gives rise to the scheme. Then by taking 
the exterior differential of Ezmdpt'^^ in the function space, it follows 

Therefore, the difference discrete symplectic structure preserving law 

At(d2W^ A Jd^^'^)) = (175) 

holds if and only if the discrete Euler-Lagrange form is closed: 

dE^^.pt'-''' = 0. (176) 

Furthermore, due to the exact forms do not change the closed condition, the null form 
may be redefined by adding certain exact forms so that the scheme may be generalized 
to a type of schemes while the difference discrete symplectic structure preserving law is 
the same. 

It is also interesting that this issue offers a way to generalize the scheme to the high 
order ones while the difference discrete symplectic structure preserving law is the same. 



Let us consider the 4-th order symplectic scheme as follows ^2 



(177) 

J^iV,((V,i^)'^'Ji/,,JV,i/)(i(z(«+i) + ^W)). 
Introduce a new "Hamiltonian" Ti. 

n = H- ^iV,HfJH,,JW,H, (178) 

then this 4th-order symplectic scheme becomes 

V^z^") = J-iV,7^(i(z("+^) + z("))). (179) 
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The discrete Euler-Lagrange 1-form associated with this case now can be introduced: 

E.uh^'''^ = + z^^yf{JAtZ^^^ - V,7^(i(2(^-+i) + (180) 

It is easy to check that these two discrete Euler-Lagrange forms differ by an exact form: 

E.^apt'^"^ - E^.J"^ = gda, (181) 

where a = {V JH,,JV ,H is a function of ^(^(^=+1) + z^^')). 

This means that they are cohomologically equivalent. In addition, this also indicates 
that the 4-th order midpoint scheme is difference discrete variational as well. 

5.2 Difference discrete variational principle and discrete coho- 
mological approach to multisymplectic algorithm 

As in the case of both Lagrangian and Hamiltonian formalism for classical field theory, 
the scenario of the above-established difference discrete variational principle and coho- 
mological approach to the difference discrete Lagrangian and Hamiltonian formalism for 
classical field theory should be directly applied to the difference discrete versions of PDEs 
with Lagrangian and/or Hamiltonian and to the numerical schemes for the multisymplec- 
tic algorithm. 



5.2.1 Difference discrete variational principle and multisymplectic algorithm 

In this part of the subsection, we first point out that the difference discrete field 
equations derived by the difference discrete variational principle themselves offer numerical 
schemes for the multisymplectic algorithm. We also derive a new scheme for the difference 
discrete Hamiltonian formalism by the difference discrete variational principle with the 
generalized modified Leibniz law ( |167| ).Then we study how to apply the difference discrete 
variational principle approach to the numerical schemes for the Hamiltonian-like PDEs 
named in [|^. A type of so-called Hamiltonian PDEs proposed first in 0] (see also [pT| ) 
are in fact a type of the Hamiltonian-like PDEs. 

As was just mentioned, the difference discrete Euler-Lagrange equations (|67D derived 
via the difference discrete variational principle in the difference discrete Lagrangian for- 
malism for classical field theory, i.e. 

offer themselves a numerical scheme that preserves the multisymplectic structures for the 
2-dimensional spacetime/space in difference discrete version. 

For the difference discrete Hamiltonian formalism of the field theory, the difference 
discrete canonical equations ([75|) derived via the difference discrete variational principle, 
i.e. 



a{i,j) _ '^riP 
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also offer themselves a set of numerical schemes for the multisymplectic algorithm. It has 
been proved in the subsection 4.2 that they preserve the multisymplectic structures as 
well for the difference discrete spacetime/space. 

In the last subsection, it was shown that by the difference discrete variational principle 
with the generalized modified Leibniz law the Euler midpoint-scheme had been derived 
from a difference discrete action functional. Similarly, the same generalized modified Leib- 
niz law ( |167|) can be employed to construct the numerical schemes for the multisymplectic 
algorithm. Let us consider the case in the difference discrete Hamiltonian formalism for 
the difference discrete field theory. 

We introduce the difference discrete conjugate momenta 

where 

is the difference discrete Lagrangian density with suitable discrete variables ■u°'(*+i'-?+i), 
Atu"^''^+^\ A^M"(*+i^'), and 

^"(i+ii) .- + := + u"^^'-'')), (183) 

4 

The difference discrete Hamiltonian density via the difference discrete Legendre transfor- 
mation 

Hd^'''^ := Hd(m"^*+^'^'+^\ 7r^^*+^'''\ A,u"(*+ij')) (184) 
= 7r^(^+5'^')AfM"(*'^'+5) - £^(M"(^+ij'+^), AtM"(*'^'+5), A^m"(*+5J)). 

Now the difference discrete action functional is given by 

Sd = J:^,jezCD'~'''\ (185) 

The variation of So can be manipulated in the framework of multiparameter differential 
approach and the differential of C^^^'^^ with respect to can be calculated to get 

dCo^''^'^ = (AiM^(*'^+5) _ (^)(M))^^^(^+ij) (186) 
where the notation 

vr;,(--^) = i(4(-+n4(-)) = (^)(-) (187) 



and the generalized modified Leibniz law ( 167 ) for the both Aj and A^. have been adopted. 

Now the stationary requirement against all variations along any direction of difference 
discrete variational principle leads to a new scheme, which may be called the midpoint 
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box scheme, in the difference discrete Hamihonian field theory, with assuming variations 
of both u^^'^'^+i) and vanish at infinity, as follows: 

= (^){.,.) (188) 



together with 



In fact, it is shown that this midpoint box scheme can be derived via the difference 
discrete variational principle with the generalized modified Leibniz law ( |167| ). 

Finally, let us derive the midpoint box scheme for a type of PDEs by means of the 
difference discrete variational principle with the generalized modified Leibniz law ( |167| ): 

MZi + eKZ, = V,^(Z), (189) 

where the same notations in HI have been used and = — K, = — L, e = ±. 



It can be shown |16] that there is an action functional for the PDEs ( |189|) 

S = ( d\C, (190) 



C = ^Z''{MZt + eKZ^)-S{Z), 

where C is the Lagrangian density, and the PDEs (|189| ) can be reached by the variational 
principle of action functional as the Euler-Lagrange equations. 

Let us now consider its difference discrete version formulation. First, introduce the 
difference discrete Lagrangian density as follows 

= /:^(Z('+5.^'+5), AtZ(*'^'+5), A,Z('+5.i)), (191) 

where Z*^*"'"^'-'^^) is taken as the coordinates in the configuration space on the note (i,j) 
in the difference discrete spacetime, AtZ^^'^~^^^ and A^.Z^'''^^'^^ the ones on the tangent 
space of the configuration space. In the difference discrete variational principle, they are 
regarded as variables to be variated. Now the action functional becomes 

Sd= E ^d^'''^ (192) 
ii,j)ezxz 

The variation of Sr, along the direction f3 in the configuration space of the system can be 
calculated in the framework of multi-parameter differential approach 

d 

It is straightforward to get the differential with respect to for each term under the 
summation 

rf£z)^''''^ = rfZ(*+i^'+^)^(MAiZ(*'^'+^) + KA,Z('+i^') - V,^(Z('+i^'+5))) (194) 
-iAi(dZ(^'^^+5)^MZ(^'^^+5)) _ ieA,(rfZ(^+i^^)^KZ(^+i^^)). 
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Here the variables Z^'^~^2J+2\ AtZ^^'-'+a) and Aj^Z'^^+a are regarded as the variational 
variables, the multi-parameter are omitted and the generalized modified Leibniz law 



(|167|) for the differences are adopted. 

Then the stationary requirement against all variations along any direction of difference 
discrete variational principle leads to the midpoint box scheme, with assuming variations 
of both and Z(*+2'-?') vanish at infinity, as follows: 

MAtZ(*'^'+^) + eKA^.Z('^+i^') = V,S{Z^'+^'^+^2)). (195) 

Thus, it is proved that the midpoint box scheme for the type of PDEs in IQ can also be 
derived via the difference discrete variational principle. 

5.2.2 Cohomological approach to multisymplectic algorithm 

Similar to the case for symplectic algorithm, it should be emphasized that in the conven- 
tional approach to the multisymplectic algorithm, it is commonly accustomed to regarding 
a numerical schemes as a set of (spacetime/space-discrete) mappings. In order to justify 
whether a set of given numerical schemes are multisymplectic, the usual approach to the 
multisymplectic algorithm is to verify this set of mappings are multisymplectic preserving 
and the verification is always carried out in the solution space of the schemes (see, for 
example, @][0]). In the difference discrete version of the Euler-Lagrange cohomological 
approach, however, instead of working on the solution space, it is first to release the 
schemes away from the solution space, even if the solution space does exist, and introduce 
some suitable difference discrete Euler-Lagrange 1-forms associated with the schemes such 
that the null difference discrete Euler-Lagrange 1-forms give rise to the schemes. Then by 
taking the exterior derivative of the difference discrete Euler-Lagrange 1-forms to investi- 
gate whether it leads to a spacetime/space-discrete multisymplectic structure preserving 
law. 

We first show the multisymplectic property of the midpoint box scheme ( |188| ) in the 
cohomological approach. 

Introducing the difference discrete Euler-Lagrange 1-forms 

EJ''^^ := {Atu''^''^+-2^ - {^f^^))dTXp^'+y\ (196) 

:= (-A,7r,(^'^-) - (^)(^'^-) + A.vr;,(^-^-))rf«/^(^+i^-+^), 

and difference discrete canonical 1-forms 

^o(.j) _ T,^ii,^)du^^'^^^-2\ e'^''^'^ := 7i'/''^du''^'+'2'^\ (197) 
the equation (|186|) can be rewritten as 

dCo'-''^^ = eJ'^'^ + eJ''^'^ + At9^^''^^ - AJ^^'^^\ (198) 

Taking the exterior differential again, due to the nilpotency of d, d'^Co^^'-'^ = 0, it follows 
that 

d{E^ + E^f'^^ + Atcu°(^'^') - A^u^^''^^ = 0, (199) 

where 

^o(^J-) := rf7r^(*J) A uj^^^^i) ■= dn'/''^ A du^^'-^'^'^\ (200) 
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Thus, the necessary and sufficient condition of the multisymplectic structure preserving 
law 

A,cj°(^'^') - A^cji(^'^') = 0, (201) 

for the midpoint box scheme is that the associated difference discrete Euler-Lagrange 
1-forms are closed: 

d{Eu + E^Y''^^ = 0. (202) 

Since the null difference discrete Euler-Lagrange 1-forms corresponding to the midpoint 
box scheme and they do satisfy the closed condition so that the midpoint box scheme is 
multisymplectic. 

Let us now consider the midpoint box scheme for the type of PDEs in 0] and show 
that it is multisymplectic in the cohomological approach. 

As before, we first define the difference discrete Euler-Lagrange 1-form and two dif- 
ference discrete canonical 1-forms 

E^J'^^^ := rfZ('+5J+5)r(MAtZ(*'^'+5) + eKA,Z(^+5J) - V,5(Z('+5.i+5))), (203) 

:= ^(rfZ(*'^'+^)^MZ(*'^'+^)), e''-''^^ := ^{dZ^'+y^^KZ^'+-2'^^) (204) 

Then the differential of Cd^^'^^ in the equation ( |194| ) can be rewritten as 

dCD^''^^ = E^J''^^ - At^°(^'^) - t/\J^^''^\ (205) 

Taking the second exterior differential of the above equation, due to the nilpotency of d, 
rf2£^(i,i) = 0, it follows that 

= dE,J'^^ - Aicu0(^'^) - e^^J''^\ (206) 

where lu^^'^'^\ fj. = 0, 1 are two symplectic 2- forms 

:= -i(rfZ(*'^'+^)^ A MrfZ(''^'+^)), (207) 

Thus, the necessary and sufficient condition of the multisymplectic structure preserving 
law 

A,u;0{i,i) _ ^^^JiJ) = 0, (208) 

for midpoint box scheme of the type of PDEs in is the associated difference discrete 
Euler-Lagrange 1-form is closed: 

dEhoJ'^^^ = 0. (209) 

Since the null difference discrete Euler-Lagrange 1-form corresponding to the box scheme 
and it does satisfy the closed condition so that the midpoint box scheme for the PDEs in 
m is multisymplectic. 

It is clear that for the both midpoint box scheme for Hamiltonian field theory and for 
the type of Hamiltonian-like PDEs the discrete multisymplectic structure preserving law 
holds in function space with the discrete closed Euler-Lagrange condition in general and 
can also be required in the solution space in each case. 

Finally, it should be mentioned that the both above difference discrete variational 
principle and cohomological scenario should also be applied to other numerical schemes 
in order to justify whether they are variational and multisymplectic. 



38 



6 Concluding remarks 



A few remarks are in order: 

1. The difference discrete variational formalism widely employed in this paper is 
different from the one of the Veselov type for the discrete classical mechanics ||24 . 



It has been emphasized that the difference as discrete derivative is an entire geometric 
object. The discrete integrants can also combine together in certain manner as a geometric 
object to construct some numerical schemes as was shown in the section 5. This is more 
obvious and natural from the viewpoint of noncommutative geometry. In the continuous 
limit, the results given here by the difference discrete variational principle lead to the 
correct continuous counterparts. 

It is shown that the difference discrete variational principle works for the difference 
discrete version of classical mechanics and field theory in the both Lagrangian and Hamil- 
tonian formalisms that present themselves as symplectic or multisymplectic numerical 
schemes and furthermore for other numerical schemes in both symplectic and multisym- 
plectic algorithms respectively. And the role-played by the different Leibniz laws for the 
differences are quite important in constructing the numerical schemes. As a matter of 
fact, all numerical schemes in symplectic and multisymlectic algorithms should be derived 
by virtue of the difference discrete variational principle together with the suitable Leibniz 



law for differences. We will publish our further investigation on this issue elsewhere p6 



2. The cohomological approach adopted in this paper to the symplectic and multi- 
symplectic geometry and their difference discrete versions in the both Lagrangian and 
Hamiltonian formalism for the classical mechanics and field theory had been missed in 
other approaches (see, for example, The Euler-Lagrange cohomolog- 
ical concepts and their difference discrete versions, such as the Euler-Lagrange 1-forms, 
the null Euler-Lagrange 1-forms, the coboundary Euler-Lagrange 1-forms as well as the 
closed Euler-Lagrange conditions and their difference discrete versions, have been intro- 
duced and they have played very crucial roles in each case to show that the symplectic 
and multisymplectic structure preserving properties. It has been show that the necessary 
and sufficient condition for symplectic and multisymplectic structure preserving property 
in each case is the related closed Euler-Lagrange condition being satisfied. Therefore, 
these symplectic and multisymlectic structure preserving properties hold in the function 
space with the relevant Euler-Lagrange condition in general rather than in the solution 
space only. Although either the Euler-Lagrange equations and the canonical equations or 
the different difference discrete versions of them do preserve the relevant symplectic and 
multisymplectic structures. 

It should pointed out that the content of the Euler-Lagrange cohomology and the roles 
played by the cohomology in each case should be further studied not only in classical level 
but also in quantum level. And needless to say, this cohomological scenario should also 
be tried to apply to other kinds of the mechanics and field theories such as the ones with 
different types of constraints and so on so forth. 

3. As was shown in the text, the variational principle/difference discrete variational 
principle and the cohomological approach form a connecting link between the preceding 
and the following in either continuous or difference discrete case. And the multi-parameter 
differential approach provides a common framework for both of them. 

It has been emphasized that both the variational principle including difference discrete 
variational principle and the Euler-Lagrange cohomological approach including its differ- 
ence discrete version can be directly applied to the ODEs and PDEs and their discrete 
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versions, which offer themselves certain numerical schemes in the symplectic and multi- 
symplectic algorithms, no matter whether there are known Lagrangian and/or Hamilto- 
nian associated with. In fact, the action functional may be constructed for certain types of 
ODEs and PDEs. Thus, the variational principle/difference discrete variational principle 
and the scenario of the cohomological approach are also available. 

As was emphasized, in the cohomological approach it is always to release the ODEs, 
PDEs and numerical schemes away from their solution spaces and to work on the relevant 
function space rather than on the solution space even if it does exist. In the standard 
or conventional approaches to the numerical schemes in symplectic and multisymplectic 
algorithms, however, in order to show whether a given scheme is symplectic or multisym- 
plectic, it is always working on the solution spaces. The implication of this difference is 
quite clear. 

4. Some simple noncommutative differential calculus on the regular lattices are em- 
ployed in our approach. Since the base space coordinates t or xs are difference discretized 
and differences do not satisfy the ordinary commutative Leibniz law for the differential, 
in order to study the symplectic and multisymplectic geometry in these difference dis- 
crete systems it is natural and meaningful to make use of the noncommutative differential 
calculus. 

5. The difference discrete version of the classical mechanics and field theory can be 
further generalized. What have been studied so far is the time discrete classical mechanics 
and spacetime/space discrete classical field theory, while the configuration spaces at each 
moment of tk in the discrete mechanics and the one at each node of the discretized 
spacetime/space on the lattice and so on are still continuous. As a matter of fact, these 
may be generalized to the case that the configuration spaces at each moment and/or node 
are also discretized. For the Hamiltonian mechanics, this is closely related to the case 
of difference discrete phase space approach to the systems with separable Hamiltonian 

mm- 

6. Finally, it should be mentioned that there are lots of other problems to be further 
investigated. 



Appendix 

We have presented the noncommutative differential calculus on the regular lattice 
and its applications to the discrete symplectic algorithms with difference discrete phase 
space |]l3|T^[0. In this appendix we briefly recall some content of the noncommutative 
differential calculus on L". General approach to the noncommutative differential geometry 



can be found in pO|. 



A.l. An noncommutative differential calculus on an Abelian discrete group 

Let G be an Abelian discrete group with a generator t, A the algebra of complex 
valued functions on G. 

The left and/or right multiplication of a generator of G on its element are commute 
to each other since G is Abelian. Let us introduce right action on A that is given by 

RJ{a) = f{a-t), (210) 

where f & A, a E G, t the generator and ■ the group multiplication. 
Let V be the space of vector fields, 

V = span{dt}, 
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where dt is the derivative with respect to the generator t given by 



{dtf){a) = RJ{a) - /(a) = f{a ■ t) - f{a). (211) 
The dual space of V , the space of 1-form, is = span{x^} that is dual to V: 

X\dt) = 1. (212) 
The whole differential algebra Q* can also be defined as fi* = with A = ^2°. 

n=0,l 

Let us define the exterior differentiation in Q* such that d : Q'^ ^ Q^. It acts on a 
0-form f & uj^ = A is as follows 

df = dJx' e n\ (213) 

Now, the following theorem can straightforwardly be proved. 
Theorem: The exterior differential operator d is nilpotent and satisfies 

(a) {df){v)=v{f),veV,feQ', 

(b) d{uj(g)uj') = duj(g)uj' + {-lY^^^'^^uj(g)duj', uj,uj'en% ^ ' 

if and only if 

(1) ® X' = {-i)x' ® x\ 

(2) dx' = 0, (215) 

(3) x'f = {Rtf)x'- 

This theorem indicates that the (g)-product should be defined as A-product and x*, the 
base of fl^, can be denoted as dt. 

As was shown here, in order to establish a well-defined differential algebra, it is nec- 
essary and sufficient to introduce the noncommutative property of the multiplication 
between function and 1-form. 

The conjugation * on the whole differential algebra ^2* and metric on discrete Abelian 
group can also be defined. 

In the case that the Abelian group is a discrete translation group with one generator 
on R^, the action of the group generates a discrete chain with equal step- length. In 
the function space on L^, the integrals can be defined (see, for example, ||25|) as follows: 

/d^oo df{t) := Ekez ^tfitk) = /(4=+oo) - /(4=-oo), (216) 



A. 2. An noncommutative differential calculus on Regular Lattice 

Let us consider the discrete translation group = ^^-^G^ with m generators, A 
the function space on G"^ and a regular lattice with equal spacing in each direction of 

on an m-dimensional space i?™. Here the i-th discrete translation group with one 
generator acting on 1-dimensional space with coordinate q in such a way: 

Rg^ ■■ q'n ^ 1n+l = 1n + h\ G (217) 

Rqi the discrete translation operation of the group and it maps of n-th node of 
g' to the one g^^^ at n + 1-th node, the discrete translation step-length along this 
direction and i?+ the positive real number. It is easy to see that the action of G* on i-th 
1-dimensional space R^ generates the i-th chain L*, -i = 1, ■ ■ ■ ,m, with equal spacing K^. 
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Similarly, the regular lattice with equal spacing /i* on each direction is generated by 
acting on R^. Since there is a one-to-one correspondence between nodes on U and 
elements of G*, one may simply regard U as G*. For the same reason, one may simply 
regard as G^. 

On the nodes of the regular lattice L™, there are discrete coordinates ql^), i = 1, ■ ■ ■ , m. 
There is a set of generators in the discrete translation group G"^ acting on L"^ in such a 
way: 

Rq^ ■ ^ Qn+i, i = l,---,m. (218) 

With respect to the generators there is a set of independent derivatives dgi on fn{(t) = 
fiqli) G A. They should be defined as the correspondent forward differences of the 
functions valued at two nearest nodes, i.e. 

= A,,/(g;) = - = - /(g;)]. (219) 

The differential 1-form is defined by 

df = d,.fdq^ = Ag.fdq\ feA. (220) 

The 2-forms and the whole differential algebra Q* can also be defined. Here d is the 
exterior differential operator. Similarly, the following theorem can be proved for d. 
Theorem.- d is nilpotent and satisfies the Leibniz rule, i.e. 

^2 = 0, d{uj Aid') = dcu Au' + {-IY^sHlu Adco', uj,uj'en*, (221) 

if and only if 

f{q^ + h^)dq' = dq'f{q'), 
q^dq^ — dq^q^ = —h^dq\ 

The above two equations show the noncommutative properties between the functions 
(including the coordinates) and differential forms. 

From these properties, it follows the modified Leibniz rule for derivatives (i.e. the 
forward differences): 

A,.(/ • g) = A,./ ■ g + {R^.f} ■ A^.g. (223) 

It should be noted that first from the definition of the forward difference, the more 
general Leibniz law may holds. For example, the generalized modified Leibniz law (|167|) 
used in the section 5 is the special case of a = 1/2, while the above modified Leibniz law 
( |223|) is corresponding to a = 1, in the following general Leibniz law with an arbitrary 
parameter a G [0, 1] 



Mf^'^ ■ 9^'^) = Hl^/^'""'^ + (1 - ci)f^'^)Atg^'^ + A,/W((l - a)(7(^+i) + a^7('=))}|a=i/2 
= -{Atf^^^ ■ ^(^=+1/2) _|_ f{k+i/2) . At^W), /C^+i/a) = i/2(/(^+i) + /W). 

(224) 

It is straightforward to prove that this general one (p24D consists with the definition for 
the forward difference. Namely, 

A,(/W.^«) :=l(/('^+^)-^('=+i)-/(^)-^(^)). 

But the price has to paid is that the corresponding Leibniz law for the exterior differential 
operator d no long holds (except for the case of a = 1). 
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Secondly, the definitions and relations given above for the noncommutative differential 
calculus on the regular lattice are at least formally very similar to the ones in the 
ordinary commutative differential calculus on i?*". The differences between the two cases 
are commutative or not. 

Similarly, the contraction between forms and differences can be defined as the same 
as the one in K^: 

< dq' A dq^, Agk >= dq'5l, 

< dq' A dq^, Ag,Agi >= 6i6i, (225) 
ixf<y =< a, Xf >, / G v4', a G Q^, 

where Xf is the Hamiltonian vector field of /. 

The Hodge * operator and the co-differential operator 

on the regular lattice can also be defined similarly as the ones on i?™ (see, for example. 



|25|). Consequently, The Laplacian on the lattice L"^ may also given by 

AL = d5L + 5Ld. (226) 

It is in fact the discrete counterpart of the Laplacian A on K^. For other objects and/or 
properties on i?™, there may have the discrete counterparts on L"^ as well. For example, 
the null-divergence equation of a form u on reads 

5a = 0. (227) 

Its counterpart on the lattice is simply 

Sl^l = 0. (228) 

This is the right forward difference version of the divergence-free equation. 

The discrete version of the integral on can also be defined (see, for example, [p5[] ). 

In the case of L^'™ G i?^'"^ with Lorentz signature, these equations become the con- 
servation law of a and its difference version of a^. This is available not only for the 
symplectic geometry and symplectic algorithms but also the multisymplectic geometry 
and multisymplectic algorithms as well. It should be emphasized that for the difference 
discrete counterparts on the lattice, they obey the noncommutative differential calculus 
on the lattice rather than the commutative differential calculus on i?™. This is the 
most important point. 
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